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ABSTRACT 

We have obtained a position (at sub-arcsecond accuracy) of the submillime- 
ter bright source GOODS 850-5 (also known as GN10) in the GOODS North 
field using the IRAM Plateau de Bure interferometer at 1.25 mm wavelengths 
(MM J123633+6214.1, flux density: Si. 25 mm = 5.0 ± 1.0 mJy). This source has 
no optical counterpart in deep ACS imaging down to a limiting magnitude of 
i 775 = 28.4 mag and its position is coincident w ith the position fo und in recent 



sub- millimeter mapping obtained at the SMA (IWang et al.l 120071 ). Using deep 
VLA imaging at 20 cm, we find a radio source (S20 cm = 32.7 ± 4.3 fiJj) at 
the same position that is significantly brighter than reported in Wang et al. The 
source is detected by Spitzer in IRAC as well as at 24 /mi. We apply different pho- 
tometric redshift estimators using measurements of the dusty, mid/far-infrared 
part of the SED and derive a redshift z ~ 4. Given our detection in the millime 



ter a nd radio we consider a significantly higher redshift (e.g., z ~ 6 IWang et al. 



20071 ) unlikely. MM J123633+6214.1 alias GOODS 850-5 nevertheless consti- 



tutes a bright representative of the high-redshift tail of the submillimeter galaxy 



1 Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and IGN (Spain). The National Radio Astronomy Observatory is a 
facility of the National Science Foundation operated under cooperative agreement by Associated Universities, 
Inc. 
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population that may contribute a significant fraction to the (sub) millimeter back- 
ground. 

Subject headings: Galaxies: formation — Galaxies: high-redshift — Galaxies: 
starburst — Galaxies: individual (MM J123633+6214.1) — Infrared: galaxies - 
Submillimeter 



1. Introduction 



Research in the last decade showed that to understand the assembly an d formation 



of m assive galaxies it is critical to study Submillimeter Gal axies (SMGs; see iBlain et al. 



20021 . for a review). Since their first detection ten years ago (ISmail et al.lll997l ). more than 



MAMBO in the submillimeter and millimeter regime (e.g., 


Hushes et al. 


1998; 


Bertoldi et al. 


2000; 


Dannerbauer et al. 


2002. 


2004 




Smail et al. 


20C 


2: 


Greve et al. 


2004; 


Wane et al. 


2004; 


Pope et al. 


20061: 


Coppin et al. 


2006; 


Bertoldi et al. 


2007 


). Due to the large beam size in 



the (sub)mm (SCUBA: 15", MAMBO: 11"), it is impossible to do a proper identification 
in other wavebands (such as the optical and near-infrared) based on bolometer data only. 
Thus, typically interferometric observations at radio wavelength with the VLA were used to 
identify the locat i on of the SMGs, yielding a radio identification rate of ~ 50 — 70% (e.g., 



Smail et al.ll2000l : iDannerbauer et al.ll2004l ; IWang et all 12004 : IPope et al.ll2006l : llvison et al. 
20071 ). depending on both the depths of the radio and the bolometer maps. The main 



drawback of this technique is, however, that the radio flux is redshift depen dent, whereas the 
sub) mm flux density at z > 1 is not (due to the inverse K correction; e. g..lBlain fc Longair 



1993|): one would thus expect the highest redshift (z > 3) sources (e.g.. iBarger et al.ll2000l ) 
to be very fai nt at radio wavelengths if they were to follow the radio-FIR correlation (e.g., 

re ra dio-identified SMG with the highest redshift lies at 
20051 ). Based on this sample, CO observat ions of radio 



Condon! Il992). We note that t 



3.6 fjChapman et al 



2003 



identified SMGs (e.g., iNeri et all 120031 : iGenzel et all 120031 : iGreve et all 120051) lead to the 



main conclusion that SMGs are massive galaxies in formation and I Chapman et al.l (120051 ) 
argue that the contribution to the star-formation rate density (SFRD) from t his population is 



missed by U V-selected surveys. If one only uses radio-identified sources (e.g. JChapman et al. 



20031 . l2005l ) one thus has to be aware of the fact that the radio selection introduces a bias 
towards the low end of the true redshift distribution. To summarize, our knowledge about 
the nature and properties of the SMG population is mainly based on radio-identified SMGs 
— a biased sample. 



One of the open questions in the investigations of SMGs is if there exists a signifi- 
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cant very high-redsh ift tail (z > 4; e.g., iDannerbauer et al.ll2002l . 12004 ; IValiante et al.l 120071 ; 
Younger et al.ll2007l ). The existence of such a tail would n ot be in contradictio n to constraints 
from measurements of the cosmic infrared background (IGispert et al.ll2000l ). Interferomet- 
ric observations in the (sub)mm hold the promise to locate the positions of bright SMGs 
to sub-arcsec accuracy. One of the brightest SCUBA sources in the GOODS North field 
is GOODS 850-5 (Wang et al. 2004; this sour ce is labeled GN 10 in Pope et al. 20050) 



with a submm flux Sg^ n m = 12.94 ±2.14 mJy ( Wang et al. 2004 Until very recently this 



bright SMG did not have a known cou nterpart which motivat ed observations with the IRAM 
Plateau de Bure Interferometer ( PdBI; iGuilloteau et al.lll992f) to obtain an accurate position 
for this source. In the meanwhile, Wang et ali ( 20071 ) reported the location of GOODS 850-5 
through SMA observations. Based on their detection in the submm (and subsequent iden- 
tification in IR ob servations) and a 2 o detection in the radio (Si. 4 chz = 18.7 ± 8.6 //Jy), 
Wang et al.l (120071 ) argue that the redshift of GOODS 850-5 is z > 4 and could even lie 
between z = 5.6 — 8.0 - with a much lower probability at z< 4.5. 

Here we report the results based on our PdBI observations. We ad opt the cosmological 
param eters fl ma uer = 0.27, fl\ = 0.73, and H = 71 km s _1 Mpc' 1 ( jSpergel et al.l 12003 



20071 ). Furthermore, we assu me a Salpeter initi al mass function from 0.1 to 100 M & . We 
use the AB magnitude scale (jOke fc Gunnlll983l ). 



2. Observations 

The final map resulting from our observations between December 2005 and August 2006 
with the Plateau de Bure millimeter interferometer in C and D configuration are based on 
6.8 hours of on-source observi ng time. The phase center (see Tabled]) was set at the nominal 
SCUBA position published in lWang et al.l (120041 ) . The source was observed in the 1 mm and 
3 mm band simultaneously. The 1 mm receiver was tuned at 240.000 GHz (corresponding 
to 1.25 mm) to obtain an accurate position. The 3 mm receiver was tuned to 98.3188 GHz, 
i.e., the frequency of the CO(2-l) tr ansition redshifted to z = 1.3448. This is the redshift 



of the optical source GOODS850-5a (jWang et al.l 120041 ) which is situated close (~ 2.4") to 



the SCUBA bolometer position. The total bandwith of 580 MHz covers a redshift space of 
Az = 0.0138 at that redshift (i.e. within the redshift uncertainty of the source, cf: z = 1.344 
in Cohen et al. 2000 and Pope et al. 2006; z = 1.345 in Wang et al. 2004; z = 1.34476 in 
Wirth et al. 2004). 



1 The distance between these two positions is 2.8", i.e. within the SCUBA bolometer position uncertainties. 
Pope et al.l (|2005l ) report S'sso Mm = 11-3 ± 1.6 mJy 
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The data were calibrated through observations of standard bandpass (3C273, J0418+380), 
phase/amplitude (1044+719, 1150+497) and flux calibrators (CRL 618, 3C273, 1044+719, 
1150+497) and reduced with the GILDAS software packages CLIC and MAP. The FWHM 
of the beam is 1.4" x 1.1" at 240 GHz and 3.4" x 3.2" at 98.3 GHz. 

We detected the dust continuum emission of GOODS 850-5 at 5 a significance at 
S1.25 mm = 5.0 ± 1.0 mJy. The derived position is RA(J2000) = 12:36:33.45, DEC(J2000) = 
62:14:08.95 (hereafter: MM J123633+6214.1) with an estimated position uncertainty of 
±0.3". The position and flux density of MM J123633+6214.1 are derived from point source 
fits of the visibilities in the UV plane and ar e consistent with t he values derived in the image 
plane. The difference between the SCUBA ( IWang et al.ll2004l ) and our millimeter interfero- 
metric position is only -0.02 " in RA and - 0.48 " in DEC, and our position is consistent with 



the SMA position from IWang et al.l (120071 ). see summary in Table [TJ In the 3 mm band no 



line was detected down to a rms of 0.8 mJy beam -1 (channel width: 40 MHz, 120 km s -1 ) on 
the position of the optical source at z = 1.3448 (and also at the position of GOODS 850-5). 
This corresponds to an upper limit (3 a) of the integrated fl ux intensity of 1.2 Jy km s _1 , 



adopting a typical CO width for a SMG of 500 km s 1 (e.g.. lGreve et al.l 120051 ) . However, 



we note that a number of CO d etections of SMGs have even lower integrated flux intensities 
than this 3 a upper limit (e.g. 



Greve et al||2005j ). 



We have used new, deep VLA radio continuum observations (Morrison et al. 2008) 
to constrain the radio flux density of GOODS 850-5. These observations reveal a highly 
significant radio source (S1.4 ghz = 32.7 + 4.3 A t Jy); i. e. a source that is significantly brighter 
than the upper limit reported in IWang et al.l (120071). Our derived flux density is however 
in agreement with the 3 a detection reported in iPope et al.l (120061 ) . The radio emission is 
unresolved with a FWHM ~ 1.7". 

In our analysis, we use the exquisite multi-wavelength dataset available for the GOODS 
North field (e.g., Alexander et al. 2003; Giavalisco et al. 2004; Dickinson et al., in prepara- 
tion) and have corrected these data products for t he offset of -0.38 " in declination between 
VLA and the GOODS imaging products (see also lPope et al.l 120061 ). 



3. Discussion 

3.1. Multi-wavelength Properties of MM J123633+6214.1 alias GOODS 850-5 



The PdBI/VLA source is undetected in the deep ACS images down to i 775 = 28.4 mag. 
However GOO DS 850-5 is detec ted at low significance in Spitzer IRAC measurements at the 
few fj, Jy level (IPope et al.l 120061 ) as summarized in Table [2j The source is also detected at 
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24 fxm (Table [2]) but undetected both at MIPS 70 jttm and 160 fxva at a 3 a level (Sm < 
2.0 mJy and Si 60 Mm < 15.0 mJy). However, due to the large beam size (18" at 70 fim 
and 40" at 160 /im) our source would be confus ed with other obje cts within the beam. We 
conclude that the IRAC counterpart presented in lPope et al.l (120061 ) is the proper counterpart 
of MM J123633+6214.1§ 

Due to the lack of a spectroscopic identification, we apply several photometric redshift 
estimators focusing on different parts of the SED in order to put constr aints on the possible 



redsh ift-range of the source. Applying the radio/ (sub) millimeter index (jCarilli fc Yunlll999 



20001 ) on the PdBI flux and the radio flux we es timate a redshif t z= 3 . 811^. Furthermore, 
we use the infrared galaxy SED templates from IChary fc Elbaa (120011 ) to fit the SED of our 
source (using the MIPS 24 /im, SCUBA, PdBI and VLA flux) and this method proposes a 
redshift of z ~ 3.3 ± 1.0, — within the errors consistent with the radio/(sub)mm spectral 
index estimates — and infrared luminosity Lm ~ 1 x 1O 13 L . We obtain the same result 
by only using the MIPS 24 /zm and the VLA flux. The use of the 850/1. 2 mm ratio as a 
proper redshift indicator is very limited (see Fig. 13 in iGreve et al.l 120041 ) . However, the 
SCUBA to PdBI flux r atio of 2.6 limits th e possible redshift range from z ~ — 4 (see plot in 
Greve et al.ll2004f ) . As IWang et al.l (120071) w e have also applied the widely used photometric 
redshift code hyperz (jBolzonella et al.ll2000l ) on the optical/near- infrared data from ACS and 
IRAC. Our r edshift estimates g ive a redshift z ~ 6, i.e. consistent with the range of values 
proposed by IWang et al.l (120071 ). However, we note that with this method only the stellar 
component of the source is fitted and the applied templates do not reflect the characteristic 
properties of infrared galaxies. 

Assuming a redshift z ~ 4, taking into accou nt the optical/ n ear-IR measurements from 
ACS/Spitzer and applying the reddening l aw by ICalzetti et al.l (120001) we estimate t he ex- 
tinction to be Ay ~ 3. For comparison, Smail et al. ( 2004 ) and ISwinbank et al. ( 2004 ) 
derived for SMG samples mean reddening of Ay ~ 1.70 and 3.0, based on optical/near-IR 
colors. 

The source is undetected both in the soft (0.5-2 k eV) and hard (2-8 keV ) bands in 



Alexander et~al1l2003h . However 



the 2 Ms Chandra observations of GOODS North (e.g. 
the lack of X-ray detection ca nnot exclude the p resence of an heavily obscured AGN in 
GOODS 850-5. Using IRAC, llvison et al.l (120041 ) proposed a diagnostic tool in order to 
distinguish between AGN and starburst as dominating energy source. Applying this tool 
on MM J123633+6214.1 would propose an AGN at redshift z ~ 2. However, we note that 



3 The IRAC an d MIPS 24 am flux es by 
merits reported bv lWang et al.l (|2007l ) 



Pope et al 



(|2006l ) are consistent within 1.5 a with the measure- 
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Lacy et al.l (120041 ) and IStern et al.l (120051 ) pointed out that starbursts at z 
the IRAC selection for z ~ 1 - 2 AGNs. 



4 contaminate 



Using either the template fitting from lChary &: Elbazl (120011 ) (resulting in an infrared lu- 
minosity Lir ~ 1 x 10 13 Lf7) and tak ing into account the conversion between infrared luminos- 
ity Ljr and SFR (IKennicuttl 1 1 998l ) ) . we derive a star formation rate SFR « 1800 M e yr~ l . 
By simply using the P dBI flux of 5.0 mJy at 1.25 mm we obtain a star formation rate 
SFR » 2200 M yr' 1 jBertoldi et al-lbooj ). This is consistent with the SFR derived from 
the radio flux for a source at z ~ 4 of SFR w 2000 M & yr~ x (see discussion in Ivison et al. 
2002). 

We conclude that the multi-wavelength data suggests that MM J123633+6314 is a 
starburst at z ~4. 



3.2. Radio Faint SMGs in GOODS North 



A total of 21 SMGs in GOODS North are securely identified (IPope et all 12009 ). In 
addition to MM J123633+6214.1 there are four SMGs that have even fainter VLA 20 cm 
fluxes (< 33 /iJy) than our source. Similar to the counterpart of MM J123633+6214.1 
none of them has an ACS counterpart and their I RAC fluxes are faint (< 20 //Jy). One of 



these SMGs is the well-studied source HDF850.1 (IHughes et al.l Il998l : iDownes et al.l Il999 



Dunlop et al.ll2004l ) which is one of the excellent candidates to lie at z ~ 4 or higher. To 
summarize, about 20% of s ecure ly ide ntified SMGs could lie at high redshift z ~ 4 or higher. 
Dannerbauer et al.l (120021 . |2004j) and lYounger et al.l (120071 ) also reported the discovery of 
reliable candidates of redshift z ~ 4 or higher in other deep fields. 



4. Conclusion 



We have used the IRAM Plateau de Bure Interferometer to obtain an accurate position 
and flux at 1.3 mm wavelengths of the SMG MM J123633+6214.1 alias GOODS 850-5 
(Si. 25 mm = 5.0 ± 1.0 mJy). This source is coincident with a source detected in deep VLA 
radio continuum observations at high significance [S\ A ghz = 32.7 ± 4.3 £tJy). The so urce 
has also bee n detected by rec ent observations obtained at the SMA (jWang et al.l 120071 ) and 
by Spitzer (IPope et al.l 120061 ) in the IRAC and 24 /zm bands. The photometric redshift 
indicators using the MIR- and FIR SED of MM J123633+6214.1 all point towards a redshift 
of z ~ 4. Based on our millimeter and radio detection, we consider a much higher redshift 
(z ~ 6, as proposed by Wang et al. 2007) unlikely. However, only future facilities with 
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wide mm-receivers planned e.g. the upgraded PdBI or the LMT will unambiguously solve 
the question of the true redshift of GOODS 850-5 (ALMA will unfortunately not be able to 
reach GOODS North). In any case, GOODS 850-5 represent one of the few solid candidates 
of a high-z (z > 3) SMG that contribute to the cosmic SFRD at early cosmic times. 



We would like to thank the IRAM staff for help provided during the observations and 
especially Arancha Castro- Carrizo and Roberto Neri for helpful support during the data 
reduction at IRAM Grenoble. Furthermore, we are grateful to Dave Frayer for 70 and 160 /im 
measurements on th e PdBI position and E manuele Daddi for instructive conversations and 
help in applying the Icharv fe Elbazl feoOlh SED templates on GOODS 850-5. Finally, we 
appreciate the comments of the anonymous referee which clarified a number of points and 
improved our manuscript. 

Facilities: IRAM: Interferometer, VLA, HST (ACS), Spitzer (IRAC, MIPS), CXO (ACIS). 
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Table 1: Positions of the SMG MM J123633+6214.1 alias GOODS 850-5 



ID 


Instrument 


R.A. 


Dccl. 


Offset 


Reference 








(J2000.0) 


(J2000.0) 


APdBI - Other 






MM J123633+6214.1 


PdBI 


12:36:33. 45±0. 04 


62:14:08. 95±0. 3 




this paper 




VLA Counterpart 


VLA 


12:36:33. 42±0. 01 


62:14:08. 70±0.1 


0.33" 


this paper 




GOODS 850-5 


SMA 


12:36:33. 45±0. 03 


62:14:08. 65±0. 2 


0.30" 


Wane et al. 


(20071 


GOODS 850-5 


SCUBA 


12:36:33. 45±0. 73 


62:14:09. 43±5.1 


0.48" 


Wane et al. 


(2004) 


GN 10 (SMMJ123633. 8+621408) 


SCUBA 


12:36:33. 8±1. 14° 


62:14:08±8.0° 


2.63" 


Pope et al. 


2005) 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, 
arcminutes, and arcseconds. 

The given position uncertainty is the search radius for counterparts of 8" applied by iPope et al.l ((2006) . 



Table 2: Fluxes of the SMG MM J123633+6214.1 alias GOODS 850-5 



Band Unit MM J123633+6214.1 Instruments Reference 
(1) (2) (3) (4) (5) 



X — Ray 


5-2 keV 


in- 17 


ergs cm z s 1 


<2.5 


Chandra 


Alexander ct al. (2003) 




X — flat/2 


-8 fceV 


lO" 16 


ergs cm -2 s — 1 


<1.4 


Chandra 


Alexander et al. (2003*) 




-B435 






mag 


>29.1 


ACS 


Giavalisco et al. (2004) 




^606 






mag 


>29.1 


ACS 


Giavalisco et al. (2004) 




»775 






mag 


>28.4 


ACS 


Pope et al. (2006) 




z 850 






mag 


>27.9 


ACS 


^nwalisco^e^al^ ^2004) 




S3. 6 t im 






/i.Iy 


1.21±0.39 


IRAC 


Pope et al. (2006) 




S4.5 fim 






/i.Iy 


1.96±0.36 


IRAC 


Pope et al. (2006) 




S5.8 fim 






/i.Iy 


2.72±0.88 


IRAC 


Pope et al. (2006) 




Ss.o i*,m 






/i.Iy 


5.11±1.14 


IRAC 


Pope et al. (2006) 




S24.0 /im 






/i.Iy 


30.7±5.4 


MIPS 


Pope et al. (2006) 




S70 /»m 






mjy 


<2.0 


MIPS 


Huvnh et al. (2007); Fraver et al. 


(2007) 


Sl60 |im 






mjy 


<15.0 


MIPS 


Huynh et al. (2007); Fraver et al. 


(2007) 


S850 /im- 


WaTic/04 




mjy 


12.94±2.14 


SCUBA 


Wang et al. (2004) 




S850 /»m- 


Pope06 




mjy 


U.3±1.6 


SCUBA 


Pope et al. (20061 




Ss70 IB 






mjy 


12. Oil. 4 


SMA 


Wang et al. (2007) 




Si. 25 mm 






mjy 


5.0±1.0 


PdBI 


this paper 




Si. 4 G/fs 






MJy 


32.7±4.3 


VLA 


this paper 





Note. — Col. (1): Band in which flux is measured. Col. (2): Units of the flux density measurements. 
Limits are 3 a. Optical/ne ar-IR mag nitudes are measured i n a 0.2" diameter aperture and are on the AB 
system l|Oke fc Gunnlll983h . see also iGiavalisco et al.l (|2004l ). IRAC fluxes are measured in a 4" diameter 
aperture. MIPS 24 /im flux is measured within a 5.7" aperture. For the MIPS 70 and 160 /im a 18.5" and 
40.0"aperture were used. (3): Measurements for the SMG MM J123633+6214.1. Col. (4): Instruments. 
Col. (5): References. 
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Fig. I.— 20" x 20" ACS 2 850 , IRAC 8.0 /mi, MIPS 24 and radio images of the field of 
MM J123633+6214.1. PdBI contours of 1.25 continuum emission of MM J123633+6214.1 
overlaid on the ACS z 850 image and start at 3a with steps of 1 mJy. 3.0" circles are drawn 
on the PdBI position in the IRAC 8.0 /im, MIPS 24 /xm and radio images. 



